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ABSTRACT 


The analysis involves time- variant heat conduction with the liquid- solid boundary 
moving radially through the axisymmetric heat- storage material as heat is transferred 
to the flawing cooling medium. Solutions were obtained by simultaneously solving the 
integrated nonlinear differential equations through numerical means. Quasi- steady radial 
heat conduction is utilized at an instant of time. Solutions for wall temperature, gas tem- 
perature, and solidification thickness along the tube for a 3 5- minute shadow portion of an 
Earth orbit are obtained for a solar Brayton heat receiver. A simplified solution neglect- 
ing the wall-temperature variation is also derived to compare with the numerical method. 



TRANSIENT SOLIDIFICATION OUTSIDE A COOLED PIPE WITH 
APPLICATION TO A SOLAR BRAYTON HEAT RECEIVER 
by Wellington W. Hu 
Lewis Research Center 

SUMMARY 

A transient heat-transfer problem of solidification outside a cooled pipe is analyzed. 
This problem involves time- variant heat conduction with the liquid- solid boundary moving 
radially through the axisymmetric heat- storage material as heat is transferred to the 
flowing cooling medium. Heat-balance equations for the solidified material, the tube 
wall, and the cooling medium are derived. These nonlinear differential equations are 
solved simultaneously by using numerical integration. Quasi- steady radial heat conduc- 
tion is utilized at each instant of time. 

The analysis is applied to the solidification of heat- storage material in the heat re- 
ceiver of a solar Brayton cycle space-power system. For this application, argon gas is 
the cooling medium and lithium fluoride is the heat- storage material. Solutions for wall 
temperature, gas temperature, and solidification thickness are presented for a 35-minute 
shadow portion of an Earth orbit. The effect of tube-entrance correction on heat transfer 
is also studied. 

A simplified analysis that neglects the variation of wall temperature and heat- 
transfer coefficient along the tube is also developed. For a solar Brayton application, 
the simplified solution and the numerical method are compared graphically. 


INTRODUCTION 

Transient heat transfer for axisymmetric solidification outside a cooled pipe is ana- 
lyzed. This time-variant heat-conduction problem considers a liquid- solid boundary 
moving radially as heat is transferred to the cooling medium flowing in the pipe. The 
analysis is applied to the solidification of heat- storage material in the heat receiver of a 
solar Brayton space-power system. The analysis is also applicable to other solar space- 
power systems, such as thermionic and thermoelectric, and to such problems as cylindri- 


cal metal casting and ice formation outside a cooled pipe. 

One-dimensional transient heat- conduction analyses with a change of phase have beer 
developed (refs. 1 and 2) for a plane geometry. Axisymmetric heat conduction with a 
phase change considering freezing on a cylinder is analyzed in reference 3. The results 
of reference 3 are applied as local solutions at several stations along a tube. The com- 
puting time required to get a single station solution by the method of reference 3 is often 
large (sometimes 30 min). Therefore, a numerical analysis is developed in this report, 
which allows a direct solution for the whole tube in reduced computer time. 

The analysis is demonstrated by obtaining solutions for a 30-tube solar Brayton heat 
receiver for a shadow period of 35 minutes. These solutions are plotted to show the vari 
ation of wall temperature, gas temperature, and frozen- material thickness and weight 
with position along the tube and time. In addition, the same calculations were performed 
by using heat-transfer coefficients that were corrected for tube entrance effects (ref. 4). 
A simplified solution neglecting the variation of wall temperature and heat-transfer coeffi 
cients along the tube is also developed for comparison with the numerical method. 


SYMBOLS 

c specific heat, Btu/(lb mass)(°R); J/(kg)(K) 

ir 

D diameter, ft; m 

h heat-transfer coefficient, Btu/(ft^)(sec)(°R); W/(m^)(K) 
hp heat of fusion, Btu/lb mass; JAg 
k conductivity, Btu/(ft)(sec)(°R); J/(m)(sec)(K) 

L length, ft; m 

N number of tubes 

q heat-transfer rate, Btu/sec; W 

R ratio of radius of frozen boundary to outside radius of inner tube 
r radius, ft; m 

T temperature, °R; K 
V velocity, ft/sec; m/sec 

W weight, lb mass; kg 

w total mass flow rate, lb mass/sec; kg/sec 
X scaled length, nondimensional 
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x axial distance along tube length, ft; m 

2 2 

l± absolute viscosity, (lb)(sec)/ft ; (N)(sec)/m 

3 3 

p density, lb mass/ft ; kg/m 

t time, sec 

if/ scaled time, nondimensional 
Subscripts: 
av average 

B freezing boundary 
b bulk 

e entrance 

F fusion 

f film 

g fluid 

s solid 

w tube wall 

x axial distance 

1 beginning of time interval 

2 end of time interval 


ANALYSIS 

Numerical Integration Method 

The thermal- energy- storage model considered herein is axisymmetric (fig. 1). The 
model considers a time-variant heat- conduction problem with the boundary between the 
molten liquid and the solid phase changing with time. The heat release from the latent 
heat of the storage material is conducted radially through the solid phase and convected 
to the heat-transfer fluid at the inner tube surface. Both ends of the tube are assumed to 
be insulated so that heat conduction is negligible. Furthermore, since the thermal con- 
ductivity of the tube wall is much higher than that of the storage material, which is a liq- 
uid, the temperature drop through the thin wall of the inner tube is assumed to be small 
enough to be neglected. The longitudinal conduction of heat is also neglected. The prop- 
erties of the storage material are assumed to be invariant with slight temperature changes. 
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The quantity of frozen material is such that its freezing front will not reach the outer tute 
wall. A cooling fluid such as argon passed through the inner tube, with its temperature 
varying along the tube. Therefore, the heat-transfer coefficient based on the film tempez 
ature of the cooling fluid varies along the tube. 

The equation for transient heat transfer in cylindrical coordinates is 


pc 


3 T 
P dr 



( 1 ) 


For any differential time interval dr, the heat balance in any differential section dx 
along the cylindrical wall can be treated as quasi- steady radial heat conduction. The 
latent heat due to the change of phase from liquid to solid plus the sensible heat due to the 
temperature drop of the solid equals the heat transfer from the wall to the cooling fluid. 
The equation for this heat balance, derived in the appendix, is as follows 


p s h F ( 2 nr B )dx 


dr 


B 
dr 


( c p) s T BP S < 2 " r B> <ix 77 - %( c p) s 21 “*7 f B rT s dr 


w 


= 2ffk g dx 


T B- T w 


(A8) 


In 


B 


w 


where r is the radius of the inner tube, neglecting tube wall thickness, r D is the 
radius at the frozen boundary, and Tg is the melting temperature of the solid. Here 
T , T g , and r^ are functions of time r and the axial distance x. Because the tem- 
perature difference for subcooling of the solid is small, the change of the sensible heat 
of the solid phase is much smaller than that of the latent heat. For a first approximation, 
the temperature distribution inside the solid is assumed to be linear. For the more pre- 
cise case, the logarithmic temperature distribution is used in the derivation (see appen- 
dix). The final integrated equation with the linear temperature distribution in the solid 
is as follows 
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p s( C p), 


2aT (2 T b - T W; ! - T W) 2 ) - | (In Rg av >[r B <R B; 2 - R B> ,)(R B 2 + j ♦ 1) 


w 


3 k„ 


- T m 9 (R r o - 1)(R, 


l w, 2' B, 2 

_ p s h F 


j 2 + 2) + T w> 1 (R b ^ x - 1)(R b> 1 + 2)J 


: [r| . 2 (ln R| 2 - i) - b| ! (in R | 2 - l)] (A15) 


where 


At = r 2 - 

The heat balance between the thermal- storage material and the gas at any time t and any 
portion x is given in the following equation 


T n - T 

q = 2vk dx B 5 = h (27rr )dx(T - T ) 

s In R b g w w g 


(2a) 


where 


T B - T w . r w h g |T T) 

1"R B ~* s W ' 6 


(2b) 


Writing equation (2b) explicitly for T w gives 


m r w V n R b ™ 
t b < l T g 


T = — 

w r h In R, 


wg 


+ 1 


(2c) 


The heat- transfer rate through the wall at any time r to the gas is equal to the rate of in- 
crease of enthalpy in the gas 
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(3a) 


4p)„ 3T 

V g g 

N Sx 


•w - V 


where both the heat-transfer coefficient h g and the temperature of gas or sink, are 
functions of time r and length x. Integrating equation (3a) with respect to x and using 
the average value of h g (T w - T g) over increments of Ax in the difference equation at 
any time t give: 


ND 


T g,x+Ax” T g,x w 


w Ax 


( c p1 


h 

L s» 


x^w,x ^g,x^ + ^g, x+Ax^w,x+Ax " ^g, x+AxJ ^) 


A rearrangement of equation (3b) gives the following equation for the gas temperature at 
the station x + Ax: 


T 

g,x+Ax 


NffD Ax 


T +• 
S’ x 2w 


w 


( c p) 


K 


/rp rp \ i rp 

x' w,x" g,x' g,x+Ax w,x+Ax 


g 


NttD Ax h 

1 w g,x+Ax 


( 4 ) 


In equation (4), the heat- transfer coefficient is a function of temperature and, therefore, 
a function of axial position x and time r. 

At the beginning of the tube, that is, at x = 0, and any time r, the inlet fluid tem- 
perature Tg(0, t) = Tg e is constant and known. The wall temperature T^, the frozen 
radius ratio Rg 2 anc * heat-transfer coefficient h g (0, r) can be calculated by using equa- 
tions (A15), (2), and the heat-transfer equation valid to the specific application. 

The calculation is started at x = 0, with a known value of h g (0, 0) to obtain h g (0, r) 
by iteration procedures. Then at constant time r the calculation along the tube axis is 
started at x = 0. Values of T w x+Ax and T g x+Ax are estimated by using the h g ob- 
tained at the beginning of the increment. Baseci on these estimated wall and fluid temper- 
atures, equations (A15), (2), (4), and the applicable heat-transfer equation are used to 
obtain the values of T w>x+Ax , T g x+Ax , h g x+Ax r , and Rg x+Ax . Final values of 
these parameters are obtained by iteration procedures until h g is obtained within a spe- 
cific limit. 

The simultaneous equations (A15), (2), and (4) are solved by using a stepwise approx- 
imation method and are programmed in FORTRAN IV for the IBM-7094 computer. 
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Simplified Solution 


As mentioned in the previous section, the sensible heat of the solid is a small quan- 
tity for the range of Rg considered and, therefore, was neglected. This leads to the 
development of a simplified solution. The tube wall temperature is assumed to be con- 
stant and equal to the melting point of the heat-storage material (i. e. , infinite conduc- 
tivity of the frozen material is assumed). Then the fluid temperature is not a function of 
time. The heat-balance equation between the frozen front and the cold sink is 

2 Jrr B Ps dx h F = 2sr w dx h g (T w - T g ) (5) 


Let Rg = rg/r w and simplify equation (5): 


dR 


R 


B 


— - = £ — (T - T ) 

dT V>s h F W g 


Integrating equation (6) from r = 0 to t gives 


? 2h_ 

Rb ’ 2 " 1 = VVh '" w 

w's F 


(T , - TJr 


( 6 ) 


( 7 ) 


where T w is constant and equal to Tg. The heat-transfer rate through the wall to the 
gas is equal to the rate of increase of enthalpy of the gas, as follows: 


w 

N 


( C p) 


dT 


g = 


g dx 


2?rr , h (T 
w g v w 


- T ) 


( 8 ) 


Assuming the heat- transfer coefficient to be constant and equal to an average value and 
integrating equation (8) with respect to x, yield 



exp 


2?rr Nh x 
w g 

. W H . 


(9a) 


or 
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(9b) 


T - T 

T = T - w g> e 
g,x A w 

exp 


Substituting equation (9a) into equation (7) gives 

2jrr ,h x 
w g 

W ( C p)s- 

Rg(x, t) - 1 = \f/e~ X 

where 


9 2h_T 

Rb<x ’ t) ' i '77V <t »' t *.p >hip i 

r w p s n F 


2j7r w Nh g x 



(10a) 


(10b) 


X = 


2jrr Nh x 
w g 

w H 


2h (Tt, - T Jr 
^ - g B g, e y 

r w^s h F 


Application 

The particular application investigated herein is the solidification of the heat- storage 
material in the heat receiver of a solar Brayton space-power system (fig. 2). During the 
sunlight period of an eclipsed orbit, solar heat melts the storage material. A paraboloi- 
dal mirror supplies solar thermal energy by focusing the Sun’s rays through an aperture 
onto the cavity wall where the heat- storage material is stored inside concentric tubes 
(fig. 3). Argon gas passes through the inner tube and extracts heat for use by the power 
conversion system. During the shadow period of 35 minutes, the gas absorbs the latent 
heat of fusion as the storage material freezes. The transient heat conduction through the 
axisymmetric body of the heat- storage material along the tube is solved by the analysis 
just presented. The reference conditions used in the calculations for a 30-tube heat re- 
ceiver for a solar Brayton system are shown in table I. Lithium fluoride is used as the 
heat-storage material. References 4 and 5 present data for gas heat-transfer coefficients 
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under moderately high heat flux conditions. For the application considered herein, the 
Reynolds number falls in the region of transition flow. The test data for the heat- transfer 
coefficient of reference 4 are scattered, but an average curve for these data can be repre- 
sented by the following equation: 


h 


g 


0. 00080 


k f (p f v b D w \ 
H ) 



(u) 


All properties of the fluid are based on the film temperature which is 



T) 

g 


The entrance correction factor for the heat-transfer coefficient in the transition flow is 
less prominent than that in the turbulent flow. For the worst condition, the entrance cor- 
rection factor in the turbulent flow region is obtained from the data of figure 10 of refer- 
ence 4 and applied to an entrance length of the tube equal to 10 times the tube inner diam- 
eter. 

Since the sensible heat of the solid-storage material is small, the use of a linear 
temperature distribution instead of a logarithmic temperature distribution results in a 
negligible error. Therefore, equation (A 15) was used for the calculations. 

Equations (A15), (2), (4), and (11) are used to solve for h g , Rg, T w , and T g . 

These equations need to be solved simultaneously and therefore a computer is used to ob- 
tain the solution. For this purpose, the tube was divided into 50 increments. An increase 
in the number of increments from 50 to 100 along the tube axis did not change the solution 
in the fourth decimal, and neither did a time interval of less than 40 seconds. 


RESULTS AND DISCUSSION FOR A SOLAR BRAYT0N APPLICATION 
Results of Numerical Integration Method With and Without Tube Entrance Correction 

Figure 4 shows the variation of ratio of frozen radius to inner tube radius along the 
length of tube at two different times. The tube diameter is 0. 087 foot (0. 026 m), and the 
length is 12. 5 feet (3. 81m) which is divided into 50 equal increments bounded by station 
numbers from 1 to 51. The frozen depth can be calculated by the following equation: 

L f = <R b - l)r, 
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At the end of a shadow period of about 35 minutes, the frozen depth at the entrance with 
the tube entrance correction is about eight times that at the exit (hot end). This suggests 
that more thermal storage material is required at the entrance. 

The ratio of frozen radius to the inner tube radius is about 1. 5 without correction, 
while it is 2. 02 with the entrance correction for the first station at the end of 35. 5 minutes 
of shadow period of the orbit. This increase of frozen radius ratio will require more 
thermal- storage material at the entrance of the tube. However, the entrance effects, as 
shown in figure 4, diminish rapidly after the third station. In order to avoid oversizing 
the receiver, the outside tube near the tube entrance could be made to follow the contour 
of the frozen- radius curve without entrance correction (fig. 4); then after about 15 min- 
utes of the shadow period, the frozen radius at the entrance section would reach the out- 
side tube. The result is to decrease the outlet gas temperature about 2° R (1. 1 K). Such 
a small temperature variation can be tolerated by the Brayton power system. 

Figure 5 shows the variation of the wall temperature along the length of the tube at 
three different times. Without the tube entrance correction, the calculated wall temper- 
ature drop at station 1 (at the cold end) is about 12° R (6. 5 K) from the beginning at time 
zero to the end of the 3 5- minute shadow period. However, when the entrance correction 
is considered, the wall temperature at station 1 at the end of 35 minutes of shadow period 
drops about 50° R (27. 8 K) from the postulated starting value at time zero. 

The effect of the entrance correction is also displayed in figure 6, which shows the 
temporal variation in wall temperature at the tube entrance and exit. Although the en- 
trance correction increases wall-temperature drop at the tube entrance from 12° to 50° R 
(6. 5 to 27. 8 K), this effect is markedly attenuated at the tube exit, and exit wall temper- 
ature varies less than 0. 5° R (0. 28 K) as a result of the entrance correction. Figure 7 
shows the weight distribution of thermal- storage material along the length of the tube for 
several shadow times. At station 1, the entrance effect increases the freezing weight per 
unit tube length to 2. 5 times that without correction. Figure 8 shows that the variation of 
total frozen weight of lithium fluoride as a function of time is less with the entrance cor- 
rection than without the correction. Since the freezing solid grows much faster in the en- 
trance region with the correction, the solid phase transfers relatively more sensible heat 
and less heat of fusion than is the case without correction for the entrance effects. At the 
same time, the gas temperature becomes higher in the later stations of the tube with the 
entrance correction so that it absorbs less heat of fusion and thus reduces the total freez- 
ing weight. Because the temperature of gas at the exit is higher with the entrance correc- 
tion, more heat is transferred. This additional amount of heat is from the sensible heat 
of cooling rather than from the heat of fusion. 



Simplified Solution 


The results of the simplified solution are shown in figures 9 to 11 plotted with those of 
the numerical integration method for comparison. Figure 9 shows that the frozen radius 
ratio Rg obtained from the simplified solution is lower at the front stations and higher at 
later stations along the tube. In this figure, the temperature of the wall is assumed to be 
constant and the heat-transfer coefficient is assumed to be an average of the heat-transfer 
coefficients obtained at time zero along the tube axis. This average value of heat-transfer 
coefficient is 0. 001581 Btu per square foot per second per °R (32. 3 W/(m^)(K)). 

Figure 10 shows that the effect of entrance correction was to increase the frozen ra- 
dius ratio nearly 30 percent at station 1 at the end of 35 minutes of shadow period. At a 
later station, such as 26, the gas had a higher temperature due to the entrance effect. 

The wall temperature was also higher. Therefore, the value of frozen radius ratio is 
reduced by the entrance correction. The maximum difference of the frozen radius ratio 
at the tube entrance is 3 percent less in the simplified solution (also shown in fig. 10). 

Figure 11 shows variation of gas temperature along the length of the tube. At the end 
of 35 minutes of shadow period, the gas temperature at the exit decreased from time zero 
only about 1° R (0. 56 K) with no entrance correction considered for the 30- tube solar heat 
receiver. This relatively small variation in outlet gas temperature is acceptable because 
the power output of the system would vary less than 1 percent as a consequence of the 
temperature change. With the tube entrance correction, the outlet gas temperature in- 
creased about 2. 8° R (1. 5 K) at time t = 0. 

The gas temperature obtained from the simplified solution is always lower than that 
obtained by using the numerical integration method (see fig. 11). The maximum difference 
of gas temperatures is only about 10° R (5. 5 K) and occurs approximately at station 20. 

The local heat-transfer coefficient used in the numerical integration method is greater 
than the average value used with the simplified solution from the tube entrance to sta- 
tion 20. Thereafter, the average value of heat- transfer coefficient exceeds the local 
value obtained by the numerical integration method. This results from the decrease in 
Reynolds number at the tube exit where the gas is hot, because of the change in gas prop- 
erties. 


Suggestions for Design of Heat Receiver 

The curve for frozen radius ratio Rg from the transient heat analysis suggests the 
use of a tapered outside tube to store more thermal storage material in the entrance por- 
tion of the tube. For the proper distribution of lithium fluoride weight at all the stations 
along the tube, the outside tube can be made with a taper ratio of about 1. 4, with rings to 
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separate lithium fluoride into different compartments, as shown in figure 12. 

The outlet gas temperature at the end of a shadow period of 35 minutes is only 1. 4° R 
(0. 7 K) below the initial tube exit value (fig. 11). For the same mass flow and tube length; 
an increase in diameter of 12 percent results in a decrease of the exit gas temperature of 
about 3. 6° R (2 K) for the 30-tube receiver. For the same length and diameter, but an in- 
creased number of tubes from 30 to 50, the heat-transfer area per unit mass flow will in- 
crease and, therefore, the exit gas temperature will increase (from 1951. 8° to 1965. 8° R 
(1084. 3 to 1092. 1 K)) about 14° R (7. 8 K) or 0. 7 percent for the same flow situation. Simi 
larly, a 3. 6-percent reduction in tube length from 12. 5 to 12. 05 feet (3. 81 to 3. 66 m) will 
decrease exit gas temperature by 5. 4° R (3. 0 K); thus, a given fractional change in tube 
length has a greater effect on exit gas temperature than does a change in number of tubes. 

CONCLUSION 

This analysis investigated the solidification outside a cooled pipe with application to 
the receiver design for a solar Brayton cycle space-power system. The thermal- storage 
material surrounding the heat receiver tubes has an axisymmetric shape. The heat- 
transfer problem in the receiver is a time- variant heat- conduction problem with the 
liquid- solid boundary moving through the thermal- storage material during the shadow 
portion of an orbit. The solution also applies in general to a problem such as metal cast- 
ing and ice formation, with variable wall temperature and heat flux along the tube. 

From the heat- receiver analysis, the ratio of frozen radius to the inner tube radius, 
the wall temperature, and the gas exit temperature are solved by the Lewis IBM 7094 com 
puter in less than 2 minutes of computing time. At the end of the 35- minute shadow peri- 
od, the gas exit temperature decreased only 1. 4° R (0. 7 K) below the exit temperature 
value at time zero. 

The frozen radius ratio along the tube axis is less than 2, and the weight distribution 
curve suggests a tapered design for the thermal- storage device. 

A fractional increase in tube length is more effective in increasing the gas exit tem- 
perature than is an increase in the number of tubes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 2, 1968, 

120-33-07-03-22. 
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APPENDIX - DERIVATION OF EQUATIONS OF HEAT BALANCE 


BETWEEN SOLID PHASE AND COOLING WALL 


The governing equation of the transient heat transfer in cylindrical coordinates is 


p s c S.ii ril 
p Bt r dr \ dr 


r 9I = k 9 / r il 


's 9r 3r \ 3r 


Integration from the inner radius to the freezing boundary results in 


p (c ) f r — dr = k /r„ — - r — 

sK p 4 J t dr s \ B 3r w 3r 

r w \ r B r w, 


The heat balance at the freezing boundary is 


o h dl>B - t 9T 

Po^TT ^ T — 


rT dr = 


B _ 0T 


r — dr + r^Tx 


substitute equations (A3) and (A4) into equation (A2) to obtain 




rT dr = k r 
s w 
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Multiply both sides of equation (A 5) by 2tt dx 


pr B h F 2, dx + P s (o p ) s 2»r B dx ^ - P s ( ft) to dx ± £* 


3TI 


, rT dr = k_2jrr dx 

s ^P's 9r Jt s w 3r 

w ir 


w 


(A6) 

For a diff erential time interval, the heat transfer across the cylindrical wall can be 
treated as quasi-steady radial heat conduction. Then 


w 


3T 

3r 


T_ - T 
B w 


r w m r B 

1 L W> 


(A7) 


Substituting equation (A7) into equation (A6) gives equation (A2): 


dr R 

p s h F (2»r B )dx _ + 


dx* /*Y 

W s Vs< 2 ’ rr B>‘ i * 1 7 - X B rT s * 

w 


= 2gk dx T - B " / (A8) 

J'A 

v 


Eliminate 2ir dx from both sides of equation (A8). Let R = r/r , and multiply by 


dr 


/R k 

B RT g dR = -J 


k_ T„ - T. 


s B w 


dr 


r~ 111 R B 
w 


Then multiply by In Rg 
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/ R 

RT s ^ 


V T B - V 


dr 


w 


Integrating and rearranging terms give 



Integrate by parts to the following term 



(A 10a) 


Substitute in equation (A9) 
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For the quasi- steady radial heat conduction, the temperature distribution inside the solid 
cylindrical frozen material is 


T s = 


c In R + T 


w 


where 


c = 


w 


In R 


B 


Then 



dR = 



R(c In R + T w )dR = 




(All) 


By the relation of equation (All), evaluate the first two terms on the right side of 
equation (A 10a) 


16 




(A 12b) 

where 


17 



■ IIIMIIB II mini 


In 




+ In R b> t ) 


is introduced in the integration to avoid the use of a power series. 

Substitute equations (A 12a) and (A 12b) into the equation (A9a) and multiply by 4. 
Then equation (A 10b) becomes 



- T 


w, 2 


r 1 R^ 

R I, 2 - 1 + 1 R B, 2> - *“<ta R B, 1>J - 4 B ; n 2 / R f 1 1 - < ln R B, 2 + R B, 1> 

' b 4 v 


R?> o - R? 


+ T w Jr 2 - - 1 - I |"ln(ln R r 9 ) - ln(ln R r 1 jl + - B > 2 , B > 1 - (In 

w ’ 1 V a > 1 2 L z “» 1 J 4 ln^Rgj 


R B, 2 + ln R B, Pi 


P s h 


s F 





(A 13) 


The sensible heat-drop term of the solid in equation (A12b) is very lengthy. This term is 
small compared with the heat-of-fusion term, so that the linear temperature distribution 
inside the solid phase can be assumed with negligible error. The linear relation is 


T = T + 
s w 


l B 


w 


Rr 


(R- 1) 


(A 14) 


Then equation (A9) becomes 
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I 


p ( c 

(2T B - j - T w ^ 2 ) “ ~ — ~ ( ln r b, av^ [ T B^ R B, 2 " R B, 1^ R B, 2 + R B, 1 


w 


T w , 2^ R B, 2 ~ DCRfl, 2 + + T w, 1^ R B, 1 " !>(%, 1 + 2 ) 

= [ R l, 2 ( ln R l - !) - R i, l( ln R B, 1 - X )] 


+ 1 ) 


(A 15) 
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TABLE I. - REFERENCE CONDITIONS AND PROPERTIES USED 


IN COMPUTER CALCULATIONS 


Argon a : 

Mass flow rate, lb mass/sec; kg/sec 

Gas temperature at inlet to receiver, °R; K 

Gas temperature at outlet from receiver, °R; K 

2 2 

Inlet pressure of gas, lb force/ft ; N/m 
Pressure drop, percent of inlet pressure 
Shadow portion of 96- min orbit time, min 

0.621; 0.28 
1457; 809 
1950; 1083 
1965. 6; 94 000 
2 
35 

Lithium fluoride 0 : 


Molecular weight 

25. 94 

Melting point, °R; K 

2020; 1122 

Density of solid at 505° R (280 K), lb mass/ft^; kg/m^ 

165; 2640 

Conductivity of solid at 2020° R (1122 K), Btu/(hr)(ft)(°R); W/(m)(K) 

5; 8.65 

Specific heat of solid at 2020° R (1122 K), Btu/(lb mass)(°R); kJ/(kg)(K) 

0.581; 2.43 

Heat of fusion, Btu/lb; kJ/kg 

450; 1047 


a Properties from ref. 7. 
^Data from ref. 8. 
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Figure 3. - Sketch of solar heat receiver. 



Figure 4. - Ratio of frozen radius to tube radius at different 
stations along tube with and without entrance correction. 
Tube length, 12.5 feet (3.81 m); tube diameter, 0.087 foot 
(0.026 m); number of tubes, 30. 
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Figure 5. - Waif temperature at different stations along tube 
axis with and without entrance correction. Tube length, 
12. 5 feet (3. 81 m); tube diameter, 0. 087 foot {0. 026 m); 
number of tubes, 30. 
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(a) Time, 0 to 160 seconds. 

Figure 8. - Variation of total frozen weight of lithium fluoride with time with and without 
tube entrance correction. Tube length, 12.5 feet (3.81 m); tube diameter, 0.087 foot 
(0.026 m); number of tubes, 30. 
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Figure 8. - Continued. 
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Figure 8. - Concluded. 
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Figure 9. - Comparison of frozen radius ratio for stations 
along tube axis for numerical integration method and 
simplified solution. Tube length, 12.5 feet (3.81 ml; 
tube diameter, 0.087 foot (0.026 ml; number of tubes, 
30; time, 35. 3 minutes. 
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Figure 10. - Comparison of frozen radius ratio with time for numerical 
integration method with and without entrance correction and simpli- 
fied solution. Tube length, 12.5 feet (3.81 m); tube diameter, 0.087 
foot (0. 026 m) ; number of tubes, 30. 
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Station number along tube length, X 

Figure 1L - Comparison of gas temperatures along tube for numerical 
integration method with and without effect of entrance correction 
on heat-transfer coefficient and simplified solution. Tube length, 
12.5 feet (3.81 m); tube diameter, 0.087 foot (0.026 m) ; number of 
tubes, 30. 
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Figure 12. - Tapered outside tube for distribution of thermal-energy-storage material, lithium fluoride. 


30 


NASA-Langley, 1968 


33 E-4539 



National Aeronautics and Space Administration 
Washington, D. C. 20546 


POSTAGE AND FEES PAIL 
NATIONAL AERONAUTICS / 
SPACE ADMINISTRATION 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 


Oo'j >01 t>y 1 3"’-s 

A { w> i l. : < C ;; F ’’ >'N 
K t !■’ T l AM> A I 1 ukf. 


a, Pi 04 GO'rfO i 
A iS r : k A 10'’ V / A F n L / 
ilAC.r, -tfv. XL X ICO 8711 < 


[ T 


L 


,*■' n \ T 


A C T i . 


T t- O . L l •' 


POSTMASTER: 


If Undeliverable ( Section 1 
Postal Manual ) Do Not Re 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



